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Abstract 

In two dimensions every solution to a nonlinear elliptic system div a{-,u, Du) = has Holder con- 
tinuous first derivatives provided that standard continuity, ellipticity and p-growth assumptions hold for 
some p > 2. We give an example showing that this result cannot be extended to elliptic systems in the 
subquadratic case, i. e. that weak solutions are not necessarily continuous if 1 < p < 2. Furthermore, we 
discuss related results for variational integrals. 

1 Introduction 

The aim of this paper is to investigate some regularity properties and the possible existence of singularities 
for vector- valued weak solutions u G PF^'P(r2,R^) of second order elliptic systems in divergence form 

— div a( - ,u, Du) — b{ - ,u, Du) in . (1-1) 

We further discuss some related results for the minimization problem of convex variational integrals 

J[w] / f{-,w,Dw)dx (1.2) 

in W^'P{il.,M.^). Here the case n,iV G N for n > 2, p G (l,oo) is considered, with f2 denoting a bounded 
domain in R". The two problems are closely connected in the following sense: provided that the integrand 
is sufficiently regular, minimizers of 5" solve the Euler-Lagrange system associated to 5": 

div DJ{ ■ , u, Du) = Dufi ■ , u, Du) in n . 

Nevertheless, exploiting the fact that the minimizer is a solution to the Euler-Lagrange system does often 
not lead to the desired results since this approach cannot distinguish between minimizers and extremals. 
Therefore, the regularity of weak solutions and of minimizers has to be discussed separately to a large 
extend. Various, by now classical results are available in the literature and helped to establish a quite 
general regularity theory for both the scalar (TV = 1) and the vectorial {N > 1) case. Furthermore, several 
counterexamples to full regularity were constructed in the vectorial case. In what follows, we give a short 
description of the known regularity theory and study its consequences, but also its limits for the two- 
dimensional case n = 2. We then demonstrate how regularity and smoothness of solutions depend on the 
integrability exponent p. 

We begin with a short overview on existing regularity results (for more details and an extensive list 
of references we recommend Mingione's invitation to the dark side f Min06] ). supposing always that the 
coefficients or the integrands are sufficiently regular and that they satisfy suitable assumptions (see (|1.3I) and 
(jl.4[) below). Since the fundamental papers of De Giorgi, Nash and Moser on solutions to single equations, 
the theory of scalar weak solutions or minimizers is by now well understood, establishing regularity in the 
sense that the gradients are locally Holder continuous, independently of the space dimension n. In the 
vectorial case instead first counterexamples of De Giorgi jDG68j and of Giusti and Miranda |GM68bj dating 
from 1968 have revealed that solutions to elliptic systems as well as minima of variational integrals may 
develop singularities for n > 3 even if the coefficients are analytic. Hence, in contrast to the scalar case, we 
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can in general expect only a partial regularity result, which means regularity outside a negligible set, which 
is called the singular set. Here regularity is always understood as (Holder) continuity of the solution (or of 
its gradient), and we introduce the set 

Reg^(?«) :— |a; G ri: w is locally continuous with Holder exponent a near 

for functions w G L^{n,'R^) and exponents a £ [0,1] (with the obvious inclusion Keg^_^{w) D Keg^^{w) 
for ai < 02). Partial regularity of the solution itself in dimensions n < p + 2 (which are referred to as 
low dimensions) is obtained via Morrey-type estimates. It is traced back to Campanato }Cam821 ICam87al 
ICam87b) and yields in particular the bound n — p on the Hausdorff dimension of the singular set. Partial 
regularity of the gradient instead was accomplished via Campanato-type estimates in general dimensions by 
various authors starting from the classical papers jMoreSi IGM68al IGM79[ IEva861 IFH85[ IGM861 IAF87| and 
resulting in optimal Holder continuity outside a set of Lebesgue measure zero. However, the counterexamples 
available in the literature still leave open the question of whether or not full regularity necessarily holds true in 
dimension n = 2. We now discuss in more detail the main ingredients for proving regularity, namely classical 
Morrey- and Campanato-decay estimates for the gradient as well as some particular features exhibited in 
the two-dimensional case. 

We are first interested in Holder continuity of weak solutions to or minimizers of (|1.2p . which 

will both be denoted by u. In view of Sobolev's embedding every function in W^i'«(rj,M^) with q > 2 is 
continuous with some (possibly small) Holder exponent. By taking advantage of the minimality resp. the 
system equation, it turns out that even if u is a priori only in W^'P{ft, R^) for some p < 2, then it indeed also 
belong to such W^'''{n, M^), provided that the a priori integrability is not too small, i. e. that p > po for some 
po G (1, 2) depending crucially on the structure constants. On the contrary, for small integrability exponents 
p e (l,Po) only the Morrey regularity theory is available, which states the equivalence Regg(M) = Reg;)j(w) 
for all A G (0, 1) and in fact guarantees that Rego(u) coincides with the whole domain ft possibly apart from 
a set of Hausdorff dimension less than 2—p (but this does not exclude singularities/discontinuities). 

In the next step a non-trivial relation between RegQ(M) and RegQ(DM) is established. We first recall the 
counterexamples jNec77| IHLN96] of Necas et al., where an integrand / is constructed which - in contrast 
to the examples mentioned before - depends only on the gradient variable, and where the solution to the 
related minimization problem (jl.2p for dimensions n > 5 (resp. its Euler-Lagrange equation for n > 3) 
is Lipschitz-continuous, but not of class C^. This example is important for two reasons: on the one hand 
this particular singular solution arises from the vectorial setting and not from an interaction effect with the 
(x, u)-dependency of the integrand or the coefficients; on the other hand it shows that in general the strict 
inclusion RegQ(£)u) C Rego(w) holds. However, we now focus on the two-dimensional situation which is 
very different: in fact, C^-regularity of solutions is well-known if the convex integrand resp. the coefficients 
of the system depend on the gradient variable, cf. Proposition 12.21 Moreover, by a simple comparison or 
perturbation argument, the regularity of the comparison solution is carried over to the solution of the original 
problem and implies Regg(M) = Rego(£'w). 

The last step is the regularity improvement for the gradient Du: the minimality property of u or the 
system equation can be used under quite general conditions (in particular for arbitrary dimension and 
arbitrary integrability exponents) to prove the equivalence RegQ(Z?M) = Reg^(Z?M) with /3 the optimal 
Holder exponent (given in terms of the regularity of the coefficients or the integrand with respect to the 
(x, w)-variables) and to further show that this regularity criterion applies ^"-almost everywhere on n. 

In conclusion, the following, straightforward strategy can be employed in two dimensions: 

Sobolcv Morrey freezing Campanato 

„ for p>pi „ , s estimates r, / \ argument / 7-1 \ estimates r, / 7-1 \ 

V. = Rego(u) = R-egAi") = Reg„{Du) = Regp{Du) 

(with /3, A £ (0,1) as above), and we now proceed to the precise statement of the full regularity results 
for minimizers of variational integrals and for weak solutions to elliptic system: Dealing with variational 
integrals we consider integrands f : il x x i?^^ — !> M subject to standard differentiability, growth and 
convexity assumptions: we require that z i— > f{-,-,z) is of class with jointly continuous second order 



Regularity vs. singularity in two dimensions 



3 



derivatives and that we have 



pN 



X R 



2N 



(1.3) 



Dzzf{x,u,z) is continuous on Q x 
iy\z\P < f{x,u,z) < L{l + \z\)P, 

z.(l + |z|)P-2|^|2 < DzJix,u,z)iX,\) < L{l + \z\r-^\X\\ 
\Dzf{x,u,z) - DJ{x,u,z)\ < L(l + \z\)P-^oJai{\x-x\ + \u-u\) , 

\f{x,U,z)-f{x,U,z)\ < L {1 + \z\)P UJaiilu - U\) , 

for all x,x e n, u,u e M^, z,X £ M^tv^ r^^^ ^j^-j^ g^ed L > ly > 0, ai,a2 e (0,1). Here ujp:R+ R+ 
denotes for arbitrary /3 G (0, 1] the modulus of continuity 

ujp{t) = min{l,t^}. 
Then the following full regularity result holds: 

Theorem 1.1: Letp G (1, oo) and suppose that 51 C is a bounded domain. There existspo — Po{N, L) < 
2 such that the following statement is true: whenever u G W^'P{Vl,M.^ ) is a minimizer to (jl.2p under 
the assumptions (|1.3p . then Du is locally Holder continuous in the interior of 51 with optimal exponent 
P :— min{ 23^, ^} < ^, i. e. u G C^^^ {il, M^), for p > po- The same assertion remains true for all p > 1 
if the integrand is independent of u, i. e. f{x,u,z) = f{x,z), or if u £ Ci^^{^,'R^) . 

For the treatment of elliptic systems we consider coefficients a: 51 x x for which we 

impose similar assumptions concerning differentiability, growth and ellipticity: we require that we have 



z h^> a{x,u,z) is of class (^^(M^^^, M^^) , 

\a{x, u, z) I + \Dza{x, u, z)|(l + |2;|) <L(l + |z|)^ 

Dza{x,u,z)X- X > iy{l + \z\Y'^\X\^ , 

\a{x,u, z) — a{x,u, z)\ < L (l + |z|)^ ^ LL>a{\x — x\ 



(1.4) 



for all x, a; G 51, w, u G R , and z, A G R , with a G (0, 1). For the inhomogeneity 6: 51 x 
we assume the controllable growth condition 

\b{x,u,z)\ < L{1 + \z\)P-^ 

for aU xen,ue R^, and z G R^^. The corresponding regularity result is then given as follows: 



p2N 



pN 



(1.5) 



Theorem 1.2: Letp G (1, oo) and suppose that 51 C R^ is a bounded domain. There exists pi = Pi{N, v, L) < 
2 such that the following statement is true: whenever u G W^'P{Vl,M.^) with p > pi is a weak solution to 
p.ip under the assumptions (|1.4p and (jl.5l) . then Du is locally Holder continuous in the interior of D, with 
optimal exponent a, i. e. u G Cl;^{n,R^). The same assertion remains true for all p > \ if the coefficients 
are independent of u, i. e. a{x,u,z) = a{x,z), or if u £ Cjqj.(51,R^). 

Remark: This regularity statement is extended easily to bounded solutions of inhomogeneous systems under 
a natural growth condition and the additional standard smallness assumption on ||u||ioo. 

The first regularity result Theorem 11.11 is a special case of |BeclO[ Theorem 1.4] (using in turn the 
arguments from [KM06, Theorem 1.7]) and is recalled here in order to draw a picture, as complete as 
possible, of the topic of regularity for two-dimensional elliptic problems. The second result seems not to be 
stated explicitly in the literature - even if all parts of the proof of Theorem 11.21 are essentially known. For 
this reason we give a proof in Section[5J For the optimality in both theorem we refer to [Gro02t Example 1.1] 
and |Phi83[ Section I]. 

We now comment on the existing literature and the reason for which the distinction concerning the u- 
dependence in the above statements is made. For the moment we shall restrict ourselves to the M-independent 
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case and concentrate on systems rather than on variational integrals (for which we can pass to the Euler- 
Lagrange system). Imposing a differentiable dependence on the a;- variable, Stara |Sta71] succeeded in showing 
the existence of higher order derivatives, ending up with a global Holder continuity result for Du. In the case 
of merely Holder continuous dependence on x, arguments of Campanato |Cam82[ Section 3] reveal that every 
solution u G W'^''''{il, M^) to (|l.ip is in fact Holder continuous independently of the value of p G (1, oo). This 
corresponds to the first step described in the strategy above. Similarly, regarding fractional Sobolev spaces 
as a generalization of the class of Holder continuous functions, one further has a fractional differentiability 
result, see |Min03b) . 

For the M-dependent case and with an a priori continuous solution u, the low dimensional theory guar- 
antees in a first step Holder continuity of u. At this point the philosophy is to pass to new coefficients 
d{x, z) = a{x, u{x), z), and one is then back in the it-independent case. 

After having given the background for the full regularity results, we now proceed to the main objective of 
this paper: we want to address the problem of whether full regularity necessarily holds for all p > 1 (which 
would be equivalent to setting po,pi = 1 in Theorem 11.11 Theorem 11.21 above) or whether there might arise 
singularities - a question which was posed by Campanato [Cam82] . 

With the strategy from above, the optimal Holder regularity of Du follows only outside of the open set 
RegQ(M), i.e. in the case p e (l,Pi) outside a negligible set of Hausdorff dimension less than 2 — p, and 
it is not clear to what extend this result can still be improved. In light of the full regularity results, the 
construction of an adequate system resp. functional that might provide an example with a singular solution 
demands one of the following features: either the integrand resp. coefficients have to be less regular in 
the variable (not satisfying (|1.3l) zt resp. (|1.4|) a). or they have to depend explicitly on the m- variable (in 
case of continuous dependence we also need a construction involving a discontinuous solution). We shall 
deal with both situations, working with the function u{x) = x/\x\ G W^'P{B^\W"-) for all p G (l,ri), which 
is discontinuous in one point and appears as a prominent example in the literature: Giusti and Miranda 
|GM68b] constructed a quadratic-type functional which is minimized by u for n > 3. Moreover, passing to 
the related Euler-Lagrange equation, it is at the same time also a weak solution to an elliptic system. Taking 
advantage of this construction we will provide in Section [3] for n = 2 and every p G (1, 2) an example of a 
functional and a system with ^""-dependence on x which are solved by u. In the next step we will rewrite 
the coefficients, ending up with: 

Theorem 1.3: Let u: 3 5 _^ r2 ^-^^^ ^^^-^ ^ ^ ^ ^^^^2). Then u G W^'PiB^R"^) n 

L°°{B,M.'^), and there exist coefficients a: M.^ x R^-'^ — > M.'^^ satisfying the assumptions (|1.4|) for some 
L > v > and every a G (0, 1) such that u is a weak solution of the homogeneous system div a{u, Du) — 
in B. 

Thcorcm ll.3l is established via Proposition l3 . 31 and Theorem l3.5[ and it further yields the strict inequality 
Pi > 1. At this point we recall that pi crucially depends on the structure constants, in particular on the ratio 
L/v. As a consequence the closer the integrability exponent p is chosen to 2, the greater this ratio needs 
to be chosen; see Remark 13.41 However, whether or not there exists a minimization problem with singular 
solution under the assumptions (jl.3p remains open. 

We mention briefly that x/\x\ is frequently considered as a function taking almost everywhere values in 
the unit-sphere S"~^ C M", i. e. as a function in the space W'^'P{B", 5"^^) for p < n. Direct computation 
shows that x/\x\ is p-harmonic in the sense that it satisfies 



for all (fi G Wq'P {B'^'' ,W^) D L°° (B" ,M"). In other words, it is stationary for the p-energy /g„ \Dw\p dx under 
the constraint |w| = 1 almost everywhere (note that this case is not contained in the above regularity theory, 
but it is covered by a serial of classical papers). In fact, even minimality holds, see |HLW98| IHonOl] . 

We close the introduction with some remarks on the notation: we write Bp{xo) := {x G R" : |a; — xo| < p} 
for the n-dimensional ball centered at xq G R" with radius p > 0. The function spaces used in this paper are 
the Holder spaces C'^^", Morrey-spaces LP'°' , the Sobolev spaces and (fractional) Sobolev spaces W^'P, 
with a, ^ G (0, 1], fc G N, cr > and p G [1, 00) (see e.g. ^Ada75t Chapter 7] for the definition and embedding 




Du ■ Dip dx 




\Du\ 
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theorems for fractional Sobolev spaces). Furthermore, we shall use two abbreviations: for a bounded set 
X e M" with positive Lebesgue-measure we denote the average of a function / e L^{X) by ^ f dx, and 
for ^ e R*^ we write V{^) = (1 + \^\^)^p~^'>^^C 

2 Review of some regularity results 

We collect some regularity results which are partially already available in the literature. We will only outline 
the proofs or give suitable references. We first observe that every weak solution u e W^1'P(0,K^) to (fLT]) 
with coefficients not depending explicitly on u (or with u being a priori Holder continuous) actually belongs 
to a fractional Sobolev space. More precisely, we have 

Proposition 2.1: Let u £ M^), p e (l,oo), be a weak solution to (|f .ip under the assumption (II. 4p 

and (|1.5|) . Furthermore, we suppose the coefficients to be independent of u, i.e. a{x,u,z) = a{x,z), or 
u e C°'J(0,R^) for some 7 > 0. Then V{Du) G W"^f{n,R'^) and Du e M^) for every 

a' < a. 

Proof (Sketch): Using difference quotients techniques we can argue similarly to Mingione [Min03bl proof 
of Proposition 3.f] and |Min03al proof of Proposition 5.2], where the superquadratic case was considered, 
in order to derive the fractional differentiability from an (uniform) estimate for finite difference quotients, 
ft should be noted that no fmther assumption (such as the continuity assumption |Min03b[ (1.8)]) on the 
inhomogeneity is needed. □ 

After this higher differentiability result (which implies higher integrability via fractional Sobolev embed- 
ding) , we come to an essential ingredient needed for the application of the comparison argument in the proof 
of Theorem 11.21 namely a priori estimates for solutions of a "frozen" problem. Following |Cam821 Section 
3] we see that these solutions admit second order derivatives. Using Gehring's lemma in order to deduce a 
higher integrability result of second order derivatives (or applying a version of Widman's hole filling trick 
[Wid71| as in [SSlOi Lemma 8.2]), we thus obtain (see also |Cam87b1 Theorem 3.1] for the superquadratic 
case): 

Proposition 2.2: Let v S W^'P{Bj^{xq),W^) be a weak solution to 

div ao{Dv) =0 m Br{xo) ^VL d M" 

with coefficients ao(-) under the assumptions (jl.4P n 1 9. Then there exists < e = e{n, N,p, L,q) such that 
for every p G (0, R] we have: 

\ViDv) - {ViDu))^^^Jdx < c(|)"7 \ViDv)\^dx 



f , ,2 / /J s min{n,2+e} f „ 

/ \V{Dv)\ dx < c(j-) / \V{Dv)\ dx, 



'Bp{xo) JBr{xo) 

and both constants depend only on n,N,p,L and v. 

This proposition uncovers a peculiarity of the two-dimensional case n = 2: the solution to the comparison 
problem has Holder continuous first derivatives. This helps us to obtain in a first step a Morrey-space 
regularity result for the weak solution u to (jl.ip (see also |Cam82( proof of Theorem l.I]), which in turn 
yields the global Holder regularity result in the interior of J7 C ffi.^: 

Proof (of Theorem ll.2p : We here follow the arguments in |KM06[ Section 9] where the related result 
for minimizers of variational functionals was proved for p >2. 

Step la: Determination of pi and preliminary regularity improvement of u. Via a standard Caccioppoli 
inequality and Gehring's Lemma we first recall the higher integrability result Du G Ll^^{Q,,W^) for some 
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'^^ dx. 



exponent q > p depending only on N,p and jj. This is exploited to determine the number pi{N ,L,v) < 2 
such that q{N,p, L,v) > 2 for all p > pi- Sobolev's embedding in turn implies u £ C[^'^(r2, M^) for some 
A = A(A^, L, v) > 0. For p e (l,Pi) and general coefBcients instead, local continuity of u is assumed. In fact, 
this is equivalent to local Holder continuity u G Cj^'^(51,M^) for some A > by the low dimensional theory, 
see |Cam821 Theorem l.I]. 

Step lb: Morrey-space regularity: Du £ L^^^~^{^l^M?^) for every /i > 0. Let B2r{xo) C ^ and define 
W^'P{Br[xo),^^) as the unique solution to divao(-Du) = in Bii{xQ), where the coefBcients 
are defined by freezing via 00(2;) := a{xo, z) (we note that existence and uniqueness follow from 

standard theory for monotone operators). Using u — v G Wq'^{Br{xo),^^) as a test-function in the weak 
formulation of the comparison Dirichlet problem, we deduce the energy estimate 

/ \V{Dv)\^ dx < c{p,L,i^) I {l + \V{Du)f 

Furthermore, taking into account the growth, ellipticity and continuity assumption in (|1.4p as well as 
Poincare's inequality, we find the comparison estimate 

c-^{N,p,u) ( \V{Du) ~-V{Dv)\^ dx < [ [ao{Du) - ao{Dv)] {Du - Dv) dx 
[af){Du) ~ a{x,u, Du)\ {Du ~ Dv) dx + / b{x,u, Du) {u — v) dx 

BrX^o) J Br{xo) 

< c{N,p,L)R"^ I {l + \V{Du)f)dx. (2.1) 

JBrtixo) 

To obtain the last line, different cases need to be distinguished: on the one hand we might be concerned with 
coefficients with explicit u-dependencies and the prerequisite p > pi or u a priori continuous. This situation 
is handled via the local (7°''*' -regularity of u according to Step la. On the other hand the coefficients might 
have no explicit u-dependency, for which the above estimate holds trivially with A = 1. Combining the 
previous two estimates with the decay estimate in Proposition 12.21 we hence end up with 

{l + \V{Du)\^)dx < c{N,p,L,iy)((^Y +R"^) [ {l + \V{Du)\^)dx 

'Bp{xo) ^ Jbr{xo) 

for all p € (0,i?]. An iteration procedure (see e.g. |Gia83[ Chapter III, Lemma 2.1]) then yields: for every 
fi > Q there exists a radius Rq = Rq{h, N,p, L, 1^,0) > (independent of the center xq of the balls) such that 

[ {l + \V{Du)\^)dx < p2-M 

JBpixa) 

for all p < Rq. To conclude the desired Morrey-space embedding for Du we observe that for large radii 
p> Rq the left-hand side is easily estimated by c{p, n, N,p, L, a, \\Du\\i^p(^g g^N-^)p'^~'^. We lastly note that 
this Morrey-type estimate is in particular a further regularity improvement of u in the sense that u is locally 
Holder continuous for any exponent less than 1. 

Step 2: Continuity of Du. We next apply the comparison estimate (|2.ip and the a priori estimate for v 
from Proposition 12. 2[ and wc find 

\V{Du)-{V{Du)) \'dx 

Bp{xo) 

< cR"^ [ {l + \V{Du)\^)dx + c(^Y^' [ \ViDv)\^dx 

Jbh{xo) JBnixo) 

< c{p,N,p,L,,.,a,\\Du\\L.iB,R-)) + p'+"] R-''-' 

for every /x > 0. Now we choose i? as a power of p such that all powers of p on the right-hand side are equal, 
i. e. R = p(2-H£)/(2+aA-He) ^ ^ Hcuce, we get 



L 



Bpixo) 



\V{Du)~{V{Du))^^^^Jdx < cp^'+-^^ 



2 , . r,^.^2+^ 
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for a constant c admitting the same dependencies as above. Therefore, the exponent at the right-hand side 
is strictly greater than the space dimension n = 2 if we choose /i G (0, |^) sufficiently small. Since the 
estimate is independent of the ball under consideration, we conclude from Campanato's characterization of 
Holder continuous functions that Du is in particular locally continuous in the interior of fi. 

Step 3: Optimal Holder regularity oj Du. Standard regularity regularity (e. g. summarized for all possible 
exponents in IMin06] in Theorem 4.4 and the following characterization of the singular set) may now be 
applied, which states that local continuity of Du is in fact equivalent to local Holder continuity of Du with 
optimal Holder exponent a. We thus get the desired regularity Du G Cj°'"(ri, M^^). □ 

Remark: In fact, also global regularity estimates can be achieved in a similar way. For this purpose one 
supposes that is a domain of class C^'" and then studies solutions in the space g + Wo'^'(n,R^) with 
g € C^'"(r2,K^). Via a standard flattening and transformation procedure the problem is first reduced to 
the model situation of the unit half-ball and zero-boundary values. Then all the results above need to be 
extended up to the boundary: for the extension of the fractional Sobolev estimates in Proposition 12.11 we 
refer to the approach in |DKM07l Proposition 5.1], for the a priori estimates of the frozen solution in order 
to conclude the Holder regularity of Du to jCamSTbl, Section 6] and [Bec09, Section 3]. A (quite technical) 
combination of the interior and the boundary estimates then yields the global result. 



3 An example for irregularity 

As already explained in the introduction, the previous regularity results still leave open the questions, 
namely whether there exist systems and variational integrals in the subquadratic case p S (1, 2) which admit 
discontinuous weak solutions. The construction of such integrands or system - depending discontinuously 
on the x-variable or depending explicitly on the solution u - shall be addressed to in this last section. Giusti 
and Miranda jGM68bj succeeded in showing that in dimensions n > 3 the function x/\x\ G W^''^{B^ jM."^) n 
L°°(B^, R^) is a minimizer to a quadratic- type functional and weak solution to a quasi- linear elliptic system, 
and they thus demonstrated that discontinuities may occur in dimensions n > 3. We now take advantage of 
their construction and show that in the two-dimensional case for any p e (1, 2) the map x/\x\ S VF^'P(i?^, M^) 
is a minimizer of functional satisfying the p- growth conditions ()1.3|) . but discontinuous in the a;- variable, and 
it is also a weak solution of a homogeneous elliptic system satisfying all assumptions in (|1.4p . 

Analogously to Giusti's and Miranda's construction |GM68b] we start by defining a bilinear form on M^^^ 

via 

for p G (1,2), all u G and indices K,X,i,j G {1,2}. In what follows we shall use the convention 
A{u){z, z) = ^ ■ jg^^ 2} Alj'i^)^i^^ fo"" z,z G R^^-^. We further introduce the abbreviations 

rr f \ ^ f \ ^ 2p z-u®u 
T4z) := Tr(z) + ^ . 

Due to the symmetry of A we immediately find the following two useful identities 

2p u (X) u 
2^ 1 + \u\ 

A{u){z, z)^ Z-Z + T^{z) T^{z) . (3.1) 
We next take g : R [0, 1] as a symmetric, smooth cut-off function satisfying Ijoj < g < ll(_i.i). We set 



2n u (X) u \ 

A(.).^z + T„(.)(l+^ ) 



nig := (p-l)-i(l + sup{|.g'(s)|+2|.g"(s)|s}) > 1 

sGR 

(the only benefit of this constant will be to compensate the effects of g occurring in the convexity condition) . 
We then define an integrand f{x, z) : R^ x R^^^ — s> R via 

f{x,z) := {g{\z\^)+mgAix/\x\){z,z))^ (3.2) 
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for all a; G \ {0} x M^^^ (and with arbitrary value for a; = 0). By definition, the integrand is bounded (for 
fixed z) and O-homogeneous in the x-variable, and it also satisfies the subquadratic growth and ellipticity 
condition: 

Proposition 3.1: The integrand /(•, •) defined in (j3.2p is smooth with respect to the variable z and satisfies 
the assumptions (|f .3^ 1 - (jf .3^ ^ with constants VjL depending only on p, nig. 

Proof: The smoothness of the integrand with respect to the gradient variable is guaranteed by construction. 
The assumption (|1.3p 9 on coercivity and boundedness is easily verified by taking into account the identity 
(|3.f p . for some constant v and L depending only on p and nrig (via the choice of g). Here we already note 
that L blows up as p 2. Hence, it only remains to prove ([TSja: for this purpose we first observe 

D,,f{x, z){X, X) = p (5(|z|2) + rUg Aix/\x\)iz, z))"^ 

{gi\zf) + mg A{x/\x\){z, z)) {g'{\zf) \\f + 2g"{\zf) {z ■ \f + lUg A(x/\x\){\, A)) 

- (2 - p) (g'dzp) z-X + nig A{x/\x\){z, X)f 

for all x G K.^ \ {0} and z,A G R^^^. The second inequality in (|f .3p q then follows immediately from p.ip . 
whereas for the first one we need to take advantage of the Cauchy-Schwarz inequality, of g{s) < for all 
s > 0, and of the definition of the constant mg to infer: 

D,Jix,z){X,X) > p{gi\z\^) + mgA{x/\x\){z,z))^ 

(.g(|zn + mg A(x/|x|)(z, z)) {{mg - g'{\z\^) - 2 g"{\z\^) \z\^) A(x/|x|)(A, A)) 

-{2-p) ml A{x/\x\){z, z) A(x/\x\){\ X) 

>p{g{\z\^)+mgA{xl\x\){z,z))'^A{x/\x\){\X) > c-\p,mg) {1 + \zf)^ \X\\ 

This completes the proof of the proposition. □ 

The regularity of / also allows us to study the Euler-Lagrange system for with integrand p.2p . which 
is given by 

div [{g{\Du\^) + mg A{x/\x\){Du, Du)) {g'{\z\^) Du + mg A{x/\x\)Du)] = . (3.3) 



Due to the convexity of / minimizers and critical points of the functional (|1.2p indeed coincide. This fact 
is now exploited to determine a discontinuous minimizer u, which in particular demonstrates that both the 
singular sets of u and of Du are not empty. Moreover, by the strict convexity u is even the unique minimizer 
with respect to its own boundary values. 

Proposition 3.2: Assume VL = B C , p G (1,2), and let u: i? — !• be given by u{x) — x/\x\. Then 
u G W^'^{B,'EP') n L°°(i?,R^), and u is the unique minimizer of the functional (|1.2p with integrand /(•, •) 
defined in p.2p among all functions in the class u + Wq'^{B,M.'^). 

Proof: We start with some preliminary observations and calculations: we note that u is smooth in \ {0} 
and \u\ is bounded by 1. Furthermore, for all a; G M'^ \ {0} and every k G {1, 2} we find 

2 J2 ^^D^u, = D,\u\'' = 0- 

i6{l,2} 

We next calculate for all x G \ {0} 

Duix) = -. — : j — with \Du\ — TviDu) — \x\^^ , 

,,,,,, 21 2{p-l) x®x 

A{x \x\)Du = ^ + -j^ ( , 

\A (2-p) \x\^ 

A{x/\x\){Du,Du) ^ \Du\^ + Tv{Duf = 2\x\-'^ . 
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From the first line we immediately obtain u £ W^'P{B,M.'^) and \Du{x)\ > 1 for all x G _B\{0}, which implies 
— 0. In order to verify that u a minimizer we start by recalling the identity div(|x|~"~^a; : 







for all X G 



\ {0}. Applying this in the two-dimensional case, we hence arrive at 

div [{g{\Du\'^) + rrig A{x/\x\){Du, Du)) ^ {g'{\Du\^) Du + mg A{x/\x\)Du)] 
d r 
dx^ - 

(2-p) 



4 E 

Ke{i,2} 



xp-P {A{x/\x\)Du) 



2 

R 



p-2 f 

= 2 2 



X 

1^ 

X 2 

|a;|P+i 2-p 



2(p - 1) + X 



„|2- 



\x\- 
-2x 



'P 2 div 



1 1 



kr J 



= 0. 



Therefore, observing that the expression in the divergence on the left-hand side of the previous equality 
belongs to W^^'-'^(i?, M^^^), we derive from the integration by parts formula that m is a weak solution to the 
Euler-Lagrange system p.Sp to ()1.2p . which means that m is a critical point. The strict convexity of / then 
yields the minimization property and the uniqueness, and this concludes the proof. □ 

In particular, the Euler-Lagrange system p.3p has a discontinuous solution (with coefficients still de- 
pending only on the independent and the gradient variable). Moreover, we may also take advantage of the 
particular structure of the integrand or the coefficients above, in the sense that the x-dependence occurs 
only in terms of a;/|a:|, i.e. of the minimizer itself. Expressing the x-dependence through the known solu- 
tion (and omitting the g'-term which anyway vanishes for u) leads to the following definition of coefficients 



a{u, z) : 



by 



a{u,z) :— (^g{\z\'^) + nig A{u){z, z)^ ^ A{u)z 



(3.4) 



for all (u, z) G X M^^^. We deduce essentially from Proposition l3 . 1 1 that these coefficients have the correct 
behavior concerning growth and ellipticity, and we further prove a continuity assumption with respect to the 
It- variable: 

Proposition 3.3: The coejjicients a(-, ■) defined in p.4p are smooth with respect to the variable z and satisfy 
the assumptions (|1.4p with constants v, L depending only on p, mg and with a — 1. 

Proof: We first observe that the coefficients are smooth in the gradient variable by definition and choice 
of the cut-off function g. We now use p < 2 and the boundedness of the bilinear form A (by a constant 
depending only on p, independently of u) to find 



|a(u,z)| < c{p,mg) (l -I- \z\Y ^ 



where the constant c{p,mg) blows up as p 2. Furthermore, the boundedness and ellipticity of Dza{u,z) 
is proved as in Proposition [33] (with the slight simplification that the second derivative of g does not appear 
any more). Thus, the assumptions (jl.4p 9 and p.4p ^ hold true, and it only remains to verify the continuity 
condition (|1.4P zt: here we note that the bilinear form A{u) is differentiable with respect to u with bounded 
derivatives. Therefore, also a(u, z) is differentiable with respect to u, and a(u, z) and Dua{u, z) are bounded 
by c(p, TOg)(l -f |2:|)P~^. Distinguishing the cases |m — uj > 1 and |u — uj < 1, we hence end up with 



|a(u, z) — a{u, z)\ < c{p, mg) min{|u — m|, 1} (l + \z 



for all u,u G 
cies. 



l"^ and all z G 



5.^^^. Thus, the assumptions in (|1.4p are satisfied with the asserted dependen- 

□ 



Remark 3.4: We emphasize that the ellipticity ratio L/iy of the coefficients a(-, •) blows up as p 2 by 
definition of the bilinear form A, and this property is indeed necessary for the construction of an elliptic 
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system with a discontinuous weak solution in view of the existence of the "critical" exponent pi in Theo- 
rem [TT2] (the higher integrability exponent q > p for Du depends only on the structure data; in particular, 
the difference q — p approaches zero as the ellipticity ratio L/i> blows up). 

It is now an easy consequence of Proposition l3.2l that there exists a discontinuous weak solution - namely 
the function x/\x\ as above - to the homogeneous system related to the coefficients given by (13.41) . Hence, we 
have an example of a system satisfying the assumptions (|1.4[) and admitting a weak solution with non-empty 
singular sets: 

Theorem 3.5: Assume p e (1,2) and let u: D -B — > be given by u{x) — x/\x\. Then u e 
W^'P{B,W') n L°°(_B,R^), and u is a weak solution of the system 

Arv a{u,Du) = in B , (3-5) 

with coefficients a(-, •) defined in (j3.4p . 

Proof: We observe that the choice u{x) — x/\x\ implies a{u,Du) = Dzf{-,Du) with / taken from p.2p . 
Thus the assertion follows from Proposition 13. 2[ where div Dzf{-, Du) = was calculated. □ 

Remark: The theorem also provides an example that the fractional differentiability Du e cannot be 

obtained in the general subquadratic case for weak solutions to elliptic systems depending also explicitly on 
the solution itself. 

Remark: The question remains open whether or not there exists a discontinuous minimizer of a variational 
integral 3^[-] with an integrand satisfying all the assumptions ()1.3|) . Instead of replacing x/\x\ in the coef- 
ficients in the Euler-Lagrange equation (13. 3p by u, we could also have argued on the level of the integrand 
defining 

fiu,z) := {gi\z\^)+mgAiu)iz,z))^ 

for all (u, z) G X K^^^ and then studying minimizers of the associated variational integral. It is then easy 
to check that the function x/\x\ is still a critical point, but due to the lack of convexity of the integrand 
with respect to the u- variable this does not necessarily imply the minimization property - and hence it does 
not lead in a straightforward way to a counterexample to full regularity for minimizers. So far it is not 
clear if minimality holds (as in the quadratic case for n > 3 for which Giusti and Miranda were able to take 
advantage of the Euler-Lagrange equation) or not. 
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